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We perform a renormalization group analysis of some important effective field theoretic models for 
deconfined spinous. We show that deconfined spinous are critical for an isotropic SU(N) Heiseuberg 
antiferromagnet, if A'^ is large enough. We argue that nouperturbatively this result should persist 
down to = 2 and provide further evidence for the so called deconfined quantum criticality scenario. 
Deconfined spinous are also shown to be critical for the case describing a transition between quantum 
spin nematic and dimerized phases. On the other hand, the deconfined quantum criticality scenario 
is shown to fail for a class of easy-plane models. For the cases where deconfined quantum criticality 
occurs, we calculate the critical exponent rj for the decay of the two-spin correlation function to 
first-order in e = 4 — d. We also note the scaling relation r) = d + 2{l~tp/u) connecting the exponent 
77 for the decay to the correlation length exponent v and the crossover exponent tp. 

PACS numbers: 75.30.Kz,64.60.Cn,71.30.-|-h, 



The most remarkable incarnation of the Landau- 
Ginzburg theory of phase transitions is the one embod- 
ied by Wilson's renormalization group (RG) [H . Accord- 
ing to this point of view, the Landau-Ginzburg theory is 
uniquely determined by the effective coupling constants 
obtained by integrating out high-energy modes. In this 
way, the large distance scaling behavior of different phys- 
ical quantities is governed by the fixed points in the space 
of coupling constants. This is the so called Landau- 
Ginzburg- Wilson (LGW) paradigm of phase transitions 
0. The LGW paradigm is known to fail in a number of 
quantum phase transitions. One prominent example is 
the transition between the Neel state to a valence bond 
solid (VBS) state in a two-dimensional Mott insulator 
Q. This transition features a quantum critical point 
(QCP), which is at odds with the LGW scenario that 
would predict a first-order phase transition. The cru- 
cial observation in this context is that both phases break 
symmetries in distinct spaces: the Neel state breaks the 
SU(2) symmetry of the Hamiltonian, while the paramag- 
netic VBS state breaks lattice symmetries. A continuous 
such order-order phase transition would not be captured 
by a LGW-like point of view 0]. 

For an SU(2) Heisenberg antiferromagnet the spinous 
Za are the elementary constituents of the spin orienta- 
tion field n. We have ria — zVaZ, a = 1,2,3, where 
•z. — (zi, Z2) and CTq are the Pauli matrices. This is the so 
called CP"'^ representation of the SU(2) spins. There is 
an inherent local gauge invariance in this representation, 
since n remains invariant when the spinon fields change 
by a local phase factor, i.e., e^^^^^'Za- Thus, it 

seems to be natural to effectively describe a Mott insula- 
tor through a gauge theory coupled to "spinon matter" . 
The gauge field here is an emergent photon: it is dy- 
namically generated as a consequence of the local gauge 
invariance of n in terms of the spinon fields. Note that 
only expectation values of gauge-invariant operators can 



be nonzero, in agreement with Elitzur's theorem [5|. The 
VBS order parameter is also a gauge-invariant expecta- 
tion value, since it is proportional to (n^ • rij), with i and 
i being nearest neighbor sites in a square lattice. The 
spinous are confined in both the Neel and VBS phases. 
Indeed, spinon deconfinement would make n fall apart 
leading to a vanishing of both spin and VBS order param- 
eters. A point in the phase diagram where the spinous 
may deconfine is at a QCP, where both order param- 
eters are supposed to vanish. To actually demonstrate 
that this happens, is not an easy task. The main ar- 
gument behind the concept of deconfined quantum 
criticality (DQC) is a topological one: spinon deconfine- 
ment occurs due to a destructive interference mechanism 
between the instantons and the Berry phase J2|- This 
mechanism was recently observed numerically [6| for the 
case of an easy-plane antiferromagnet. However, in this 
case the instanton cancellation mechanism leads actually 
to a weak first-order phase transition In Ref. 01 a 
first-order phase transition in an easy-plane model of vo- 
tex loops was also found, but there Berry phase effects 
were not considered. 

One of the most important aspects of DQC is the large 
value of the critical exponent ry for the decay of the 
correlation function Q{x) = {n{x) ■ n(0)) as compared 
with the value obtained through the LGW approach. 
This correlation function is highly relevant experimen- 
tally. Therefore, it is important to be able to calculate 77 
in a systematic way. The exponent rj has been calculated 
in Monte Carlo simulations for two-dimensional Heisen- 
berg antiferromagnets with instanton suppression ^] and 
with four-spin interactions The obtained results are 
T] fa 0.7 and r] ~ 0.26, respectively. 

One of the main results of the present paper will be 
the calculation of 77 in first-order in e = 4 — d, where 
d is the dimension of space-time. This will be done for 
two different DQC regimes: (i) the Neel- VBS transition 



2 



and (ii) the phase transition between quantum spin 
nematic and dimerized phases [lo| . 

It was argued in Ref. 0] that for certain isotropic 
SU{N) symmetric Heisenberg antiferromagnets the 
quantum critical point is governed by the euclidean La- 
grangian 



N 



N / N \ 



(1) 



where the parameter is tuned in such a way that the 
system is at its critical point. The above Lagrangian 
corresponds to an Abelian Higgs model in euclidean space 
with an 0{2N) global symmetry. It can also be thought 
of as the free energy of a Ginzburg-Landau (GL) model 
with N complex order parameter fields. In this case the 
upper critical dimension is four. Thus, the RG analysis 
should be made in d = 4 — e dimensions. We define 
dimensionless couplings g — ^""^Ur and / = /i~'e^, where 
Ur and Cr are the renormalized counterparts of uq and 
eo, respectively. The RG (3 functions /3/ = fxdf /d^ and 
Pg = fidg/dfj, are straightforwardly obtained employing 
standard techniques [11|: 



Pa 



-eg - 6fg + {N + 4)3^ + 6f 



(2) 



(3) 



The above RG equations are well known in the context 
of the GL model [l^. An infrared stable fixed point 
with / 7^ is found only for N large enough, namely, 
N > 182.9. Unless TV is greater than this value, no 
second-order phase transition is predicted by this RG 
analysis. In the past this result led to the conclusion 
[1^ that thermal fluctuations turn the phase transition 
in a superconductor into a first-order one, since there the 
actual number of components is TV = 1. It did not take 
too long to realize that this result is incorrect l^, 14 1. 
Actually the large, but finite, iV result reflects the strong- 
coupling features of the TV = 1 theory, which cannot be 
captured by the RG analysis in d = 4 — e dimensions. 
This does not mean that the first-order transition cannot 
occur. It turns out that the complete phase diagram fea- 
tures a tricritical point [3] at a value of the Ginzburg 
parameter k — ujsple smaller than the mean- field GL 
value separating the type I from the type II regimes, i.e., 
K = l/y/2. Earlier calculations based on duality argu- 
ments give the value K( « 0.8/-\/2, a result rece ntly 
confirmed by large scale Monte Carlo simulations [l5|. 
The weak-coupling regime at low iV captured by the RG 
functions ^ and ^ corresponds to the one in which 
K <^ Kt- An RG analysis in fixed dimension d = 3 [l^ . 



though less well controlled than the one near four di- 
mensions, indicates that the critical value of can be 
drastically reduced if a resummed higher order calcula- 
tion is performed 17[. Other interesting effects arise as 



the number of components gets large enough for the case 
where the gauge field is compact. For example, there is 
a recent numerical evidence for a Coulomb-like phase in 
a compact abelian gauge theory coupled to multiflavor 
CP^ fields [3]. 

We can make an analysis of deconfined quantum crit- 
icality that parallels the case of superconductors in the 
neighborhood of the critical temperature. In contrast to 
the superconductor, in this case the physical value of the 
parameter N is given by A'^ = 2. Hence, we have two 
classes of (meron) vortices 0, [3| ■ The above discussion 
in the context of superconductors shows that in principle 
we can also have deconfined quantum trier iticality. 

The correlation function Q{x) = (n(a;) • n(0)) in the 
(jpN-i representation reads 



g{x) = 2(z*(a;) •z(0) z(a;) •z*(0)) 
- |(|z(a:)nz(0)|2). 



(4) 



At the deconfined QCP this correlation function scales 
as Q{x) ^ [l^]. In order to compute 77, let 

us analyse the scaling behavior of the two four-spinon 
correlation functions in Eq. (|4|. 

The scaling behavior of (|z(a:;)p |z(0)p) is obtained by 
considering the scaling dimension of the operator |z(2:)p. 
This is easily obtained by performing derivatives with re- 
spect to tq of the logarithm of the functional integral and 
doing dimensional analysis. The result is a scaling be- 
havior of the form (|z(a;)|2|z(0)|2) ~ 1 / \x\'^-'^+'^\ where 
7]^ = d -\- 2{\ — 1 / v) , with V being the correlation length 
exponent. This leads to a vanishing of those correlations 
in momentum space as p — > 0, except for the mean- field 
case where rj^ — d — 2. Indeed, beyond mean-field theory 
we have v > 2/d and thus it is clear that 774 > d — 2 
when the fluctuations are included. This result is impor- 
tant because it legitimates the softening of the CP^~^ 
constraint |zp = 1. In the critical regime we can simply 
neglect the second term in Eq. (|4]). 

Let us consider now the scaling behavior of (z*(x) • 
z(0) z(a;) • z*(0)). This correlation function is associ- 
ated with a mass anisotropy term, which is obviously 
not generated by quantum fluctuations in a SU(N) the- 
ory like the one in Eq. ([1]). However, we can con- 
sider it as a source term and compute the so called 
crossover exponent tp [2l|. The exponent 77 is then ob- 
tained by replacing \/v in the expression for 7^4 by ^p/v. 
Therefore, (z*(a:) •z(0) z(a;) •z*(0)) - l/\x\'^-'^+'^ , where 
rj = d + 2(1 — (fi/i^)- This result is obtained as follows. 

The anomalous dimensions of all quadratic op- 
erators, leading to mass anisotropy or not, can 

(2) 

be derived from a "matrix exponent" 7/^^ = 



3 



lim^^oM51n[4|^,(Z-V2)^^(Z-i/2)^,]/a^^ (here a 
summation over repeated greek indices is implied), 



where Z 



(2) 



is the renormahzation associated to the 



insertion of a quadratic operator and Z^ji the spinon 
wave function renormahzation ^21]. From the eigen- 
values of this matrix exponent we can determine both 
V and Lp or, equivalenty, 774 and 77. These eigenvalues 
are the anomalous dimensions of the composite oper- 
ators |zp and z*Z/3, with a ^ (i. We will call these 
anomalous dimensions 772 and 772, respectively. Now 
it is straightforward to use dimensional analysis to 
obtain that 774 = c? — 2 — 27/2 — d + 2(1 — \/v) and 
77 = d - 2 - 277^ = d + 2(1 - ip/v). 

For the DQC regime described by Eq. ([1]) we have at 
one-loop order 

^0,0,^6 ^ -NPap^-^sg* + (3/* - g*)Iaf3,jS, (5) 
where la/S^-yS = iSa^Sf3S + 6asSf3y)/2, Pafj,jS = Sal3Sjs/N, 

and (7* and /* are the infrared stable fixed points asso- 
ciated with the /3 functions ^ and ([2]). The eigenvalue 
772 corresponding to the eigenvector d-ys determines the 
critical exponent v as l/v = 2 + ri2- The second eigen- 
value, 772 determines ip through 772 — ip/v — 2. Explicitly, 
we have 772 = —{N + l)g^, + 3/, and 772 — — -I- 3/,. 
Therefore, we obtain to order e and for N > 182.9 the 
result 



77 



1 



18 

'n 



2.9* 



(6) 



where = {IS + N + y/N"^ - 1807V - 540)e/[2A^(A^-f 4)]. 
Using the lowest value of N for which the stable fixed 
point exists (i.e., N = 183), we obtain after setting e = 1 
the resuh 77 = 609/671 « 0.9076. The LGW result to 
order e would give rj = and a small correction to order 
e^. Note that the local gauge invariance is essential in 
order to get a value smaller than one for 77. If we consider 
the model without any gauge coupling, we obtain 77 ~ 
2 — [1 — 2/{N + 4)]e, which is an expression valid for all 
values of N, since in this case the perturbative fixed point 
exists even for = 1. This leads for N — 183 and e = 1 
to the result 77 = 189/187 > 1. 

Now we consider the scaling behavior of the spin 5 = 1 
Hamiltonian 0, [ll 



(7) 



where both nearest-neighbor couplings J and K are pos- 
itive. This Hamiltonian describes the phase transition 
between a quantum spin nematic phase and a dimerized 
phase. Recent numerical results indicate that this model 
exhibits a second-order phase transition if the ratio K/ J 
is large enough. As pointed out in Ref. [lo|, the LGW 
paradigm would in this case predict a first-order phase 
transition, at odds with the results observed numerically 



[2 2( 1 . Here we will show using the RG that in this model a 
second-order phase transition occurs for large enough N. 
The field theory of the above model was derived recently 
[lo| and is given by the Lagrangian 
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C = -ie^,xdM' + \{d^-ieoA^)-D\'+ro\n\ 



2' " 

Uq + Vq 



(8) 



where vq > and D is a complex vector with three com- 
ponents. 

To see how the second-order transition emerges, let us 
write Di — (ipi + i'ipi)/\/2, with i = 1,2,3. The local 
interaction between the scalar fields become 



int 



'^0, 2 



2n2 



(9) 



The above equation features an interaction reminiscent 
of certain classical models for frustated magnetism [j^ . 
In order to perform the RG analysis, we will consider 
a generalization of the model such that ip and ■0 have 
each N components, with the physically relevant case 
corresponding to = 3. The f3 function for the gauge 
coupling is given once more by Eq. ([2|). By introduc- 
ing dimensionless couplings g = fi~'^u and h = fi~'^v we 
obtain the one-loop /? functions: 

Pg = -eg - 6fg + (N + A)g^ + 2h' - 2gh + 6f, (10) 



-eh - 6fh -{N + 2)h^ + 65/1. 



(11) 



It is useful to analyse first the case where / = 0. The 
physically meaningful case corresponds to fixed points 
where /i > 0. The relevant fixed point in this case has 
coordinates = 2e/{N'^ + 8) and K ^ [2 - N)e/{N^ + 
8). This fixed point is stable only for A^ = 3, but then 
we would have /i* < 0, which is incompatible with the 
physical constraints of the model ^1^. 

Remarkably, for /* — ie/N a stable fixed point is found 
for A^ > 232.98: 

_ 360 - 12Ar + 2Ar2 + AT''' + (Af + 2)%/A 

■9* — ^TKJICA , OAT . OAT9. . AT?.\ ^' (1^) 



2N{6A + 8N + 8Af2 + N^) 



_ 3VA - 36 - 104A' - 21 Ar2 _ ^3 
* ~ Af(64 + 8N + 8N^ + N^) ^' 



(13) 



where A = A^^ _ 224A^3 _ 2072Af2 „ 4608A' - 22896. 
Once more, just like in the isotropic case, we interpret 
the existence of stable fixed points at large values of A^ 
as a strong evidence of deconfined quantum criticality. 
The exponent 77 is calculated similarly as before, giving 
the result 77 « 0.927 for A^ = 233 and e = 1. 

Next, we consider an easy-plane version of the model 
P]). This amounts to adding an interaction term of the 
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form t!o(|zip — |22p)^/2. Previous results on the easy- 
plane model 0, Q indicated that a second-order phase 
transition would occur. This conclusion was based on 
the analysis of the deep easy-plane limit of the model. 
This regime corresponds to a large wq such that j^il^ ~ 
|z2p- However, recent Monte Carlo simulations [y, 0| 
performed in this regime showed that the transition is 
actually (weakly) first-order. 

In order to facilitate the RG calculations it is conve- 
nient to write the complete local interaction between the 
spinons in the following form: 

Ant = y(kl|' + k2n+t«0klPk2p, (14) 

where uq — u^^+vq and wq — uq~vq. Let us introduce the 
renormalized dimensionless couplings g = u^~^ and h = 
wjjT'^ , where u and w are the renormalized counterparts 
of Mo and wq, respectively. In order to have the same 
total number of complex components as before, we will 
consider N/2 components of zi and Z2, with N even. By 
this we mean a rewriting of the interaction, such that 
the system has a 0{N) x 0{N) symmetry. The one-loop 
P function for the gauge coupling is the same as before. 
The other /3 functions are 

/3, = -eg - 6.9/ + ^g^ + + 6/^ (15) 



exist for large enough N , but none of them are stable. 
Therefore, there is no deconfined quantum criticality as- 
sociated to the model (HH). 

Summarizing, we have considered three models for de- 
confined spinons. From the three models considered, only 
the one associated with an easy-plane antiferromagnet 
does not exhibit any second-order phase transition, in 
agreement with the numerical results of Refs. and 0] ■ 
Deconfined spinons were shown to govern a second-order 
phase transition for both the isotropic SU{N) antifer- 
romagnet and quantum spin nematic systems. In both 
cases we have computed the critical exponent rj using 
the scaling relation rj = d + 2(1 — 'fi/'^) in terms of the 
crossover exponent (p and the correlation length exponent 
v. Knowledge of the exponent (p is very important in the 
study of crossover behavior and stability of frustrated 
systems. 
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f3h = ~eh~6hf+{N + 2)gh + 2h^ + 6f. (16) 

It is instructive to consider first the model for vanishing 
gauge coupling (/ = 0). In this case besides the Gaussian 
{g^, = =0) and Heisenberg [g^ — 2e/{N + S) and /i, = 
0] fixed points, we have the fixed points gi — hi — e/{N+ 
4), and (^2, /12) with g2 = Ne/{N^ + 8) and /i2 = (4 - 
N)e/ (N"^ + 8). From these only the fixed point (32, is 
infrared stable, provided iV = 3. Note that for iV = 4 we 
have the realization of the deep easy-plane limit, since 
the effective interaction multiplying |2:ip|2:2p vanishes, 
although the bare coupling wq 7^ 0. However, such a 
case does not corresponds to a stable fixed point. Note 
that for = 2 the fixed point {g2, ^12) is 0(4) symmetric 
and coincides with (171, /ii). 

Note that the gauge field fluctuations, which in this 
problem are essential, generate a |zip|2;2p term. Thus, 
we have to keep wq and look in the RG treatment 
for the stability of fixed points with h = 0. It turns out 
that for all values of N, no fixed points with h^, — and 
/* = 3e/N are found. Therefore, no second-order phase 
transition takes place in this case. This is a very signif- 
icant result, since as we have discussed, the existence of 
a critical value of N above which the transition becomes 
second-order reflects the actual behavior at lower values 
of N. The complete absence of fixed points for all N 
provides a solid theoretical explanation for the numeri- 
cal results of Refs. 7] and Fixed points with ft,* 7^ 
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